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Abstract 
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We propose a technique, using interferometry of Bose-Einstein condensed alkali atoms, for the 

detection of sub-micron-range forces. It may extend present searches at 1 micron by 6 to 9 orders 
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of magnitude, deep into the theoretically interesting regime of 1000 times gravity. We give several 



examples of both four-dimensional particles (moduli), as well as higher-dimensional particles 



. vectors and scalars in a large bulk- that could mediate forces accessible by this technique. 
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I. INTRODUCTION 



Some recent theoretical ideas point to the possibility of new physics, related to gravity, at 
the sub-millimeter regime. One is the preponderance of light gravitationally coupled moduli 
suggested by string theory 

HQ 

. Another is the possibility of large sub-mm-size dimensions 
and the particles residing inside its bulk ^, U, ^| . Yet another is suggested by the magnitude 
of the vacuum energy p. 

These ideas motivated some heroic experiments that have, in the last 

seven years, extended the search for such forces from mm down to ~ 20 microns. These 
experiments involve measuring the force between two macroscopic but small objects. A 
fundamental obstacle in searching at much smaller distances is that the size of these objects 
must be reduced and therefore the expected signal force decreases; at the same time, the 
electrostatic background Van-der-Waals force increases. 

In this paper we suggest a possible way around this obstacle by considering the interaction 
of a macroscopic system with a pure quantum mechanical system consisting of a Bose- 
Einstein condensate. The latter has a significant advantage relative to a macroscopic system: 
its de Broglie phase can be measured very precisely. In addition, it can be well controlled 
and manipulated and its electromagnetic interaction with its environment is well understood 
- both theoretically and experimentally. Because of these advantages, the technique that 
we propose may extend current bounds at 1 micron by 6 to 9 orders of magnitude, and be 
sensitive to forces as small as 1000 times gravity. The approach we describe thus explores 
a region of parameter space that is complementary to the super-micron reach of upcoming 
micro-cantilever and torsion balance experiments. 

In section 2 we update the analysis of macroscopic forces below 10 microns in theories 
with light moduli. In section 3 we consider new forces from bulk gauge fields or scalars in 
large extra dimensions, taking baryon number as an example. In section 4 we propose our 
experimental technique and estimate some of the important backgrounds. We conclude with 
section 5. 
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II. FORCES FROM LIGHT MODULI 



In string theory the parameters of the standard model depend on fields, called moduli, 
whose values determine the geometry of the extra dimensions. Moduli couple with grav- 
itational strength and typically remain massless until supersymmetry is broken. So, they 
get a mass proportional to ~ F/Mpl, where F is the scale where supersymmetry breaking 
originates. In theories with gravity-mediated supersymmetry breaking, F is (10 8 TeV) 2 and 
he moduli have microscopic Compton wavelengths. However, as pointed out in reference 
lj], in theories of gauge-mediated supersymmetry breaking, F can be as small as (10 TeV) 2 
and moduli can have macroscopic Compton wavelengths and mediate macroscopic forces of 
gravitational strength. The range and magnitude of these forces, for a variety of moduli, 
were first estimated in reference jjj for y/F in the range of 10 TeV to 100 TeV. At that 
time it seemed pointless to consider larger values of F, since they lead to moduli Compton 
wavelengths which were thought to be inaccessible to macroscopic-force experiments. In this 
paper we extend the scale of y/F up to 2000 TeV, which in turn considerably extends the 
predicted parameter space for moduli- dependent forces. 

The upper limit for the value of y/F comes from cosmology: In gauge-mediated supersym- 
metry breaking, the gravitino is the lightest supersymmetric particle, with mass ~ F/Mw,. 
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Although light, its mass still must not exceed 1 keV to avoid over-closing the universe 
This in turn provides an upper limit of 2000 TeV on y/F. 

We focus on the three classes of moduli studied in reference [1] which couple directly to 
ordinary matter: the dilaton, the gauge moduli and the Yukawa moduli. Moduli-dependent 
forces can occupy a substantially larger region of parameter space than previously indicated. 
Although much recent experimental progress has been made in the search for new sub- 
millimeter forces there is ample potentially interesting parameter space awaiting 
further exploration. 



A. Gluon modulus 

We consider here a field that couples only to the standard model gluons. The effective 
coupling is given by [j| 

= \l± G * (1) 
8tt 2 M V ' 
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FIG. 1: Experimental bounds and theoretical expectations on new forces from potentials of the form 
V(r) = — Gn mi ™ 2 (l+ae~ r//A ) below 1 cm. The projected reach of the first-round BEC experiments 
is shown as a solid green line. The solid-red line indicates the reach with an improved sensitivity 



of 10 'Hz. Experimental data are from references 
expectations are discussed in the text. 
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The shown theoretical 



where X g is an undetermined coupling constant, M is expected to be of the order of the 
string scale 5 x 10 17 GeV, and the suppression factor 87T 2 accounts for the gauge coupling 
depending on moduli only at higher-order. In this way, the coupling strength is weaker than 
that of the dilaton discussed in the following subsection. 

Considering the contribution to its mass coming from the interaction in eq. (JTJ), the 
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FIG. 2: Experimental bounds on new forces from potentials of the form V(r) = —Gn " 111712 (1 + 



ae r//A ) below 1 micron. The projected reach of the first-round BEC experiments is shown as a 
solid green line. The solid-red line indicates the reach with an improved sensitivity of lCP 7 Hz. 



Experimental data are adapted from a figure in reference 
Compton wavelength of the field is []] 



Q 



(2) 



Here F is the fermion-scalar messenger mass-squared splitting, A; is a loop-integral factor of 
order 1, and N is the number of messenger multiplets. iV < 4 is required so that all gauge 
couplings remain perturbative below the GUT scale, and a bound of \fF > 30 TeV/\/N can 
be imposed from constraints on the right-handed selectron mass and a consistency condition 
that the messengers have non- negative mass-squared 
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The coupling of to the nucleon N can be expressed by 



M PL 

where is given by 



C = Q^^-Hn^n , (3) 



\ g Mpl (N\G* u G a ^\N) _ pi /5x 10 17 GeV\ 



- 8tt 2 M m N ~ " 6A « { M J ' (4) 

The interaction in eq. (jHJ) yields a potential for exchange between particles of mass mi 
and m 2 at a distance r: 

(^tt) A; 2 + m^ 

When added to gravity, eq. (0) describes an additional attractive force: 

F(r) = -G^(l + a^) , (6) 

where a is given by Q^/Att. In reference a range of A^ 9 and a was obtained by taking 
y/kN = 1 and by varying \/F between 30 and 100 TeV and A ff 1 x M/(5 x 10 17 GeV) between 
10~ 2 and 10 2 . We here expand the range for the SUSY scale y/~F up to the limit of 2000 
TeV imposed by the gravitino problem and show the results in Figure 



B. Dilaton 



The dilaton couples to nucleons with a strength of about 80 times gravity 




leading to an inter-nucleon force of about 6400 times gravity. Since it couples to all fields 
in the theory, it is expected to receive a mass ~ F/M from its strong coupling to the 
primordial supersymmetry-breaking sector. This would make its Compton wavelength less 
than 10 -2 yum (100 TeV) 2 /F, which is too short to be experimentally observed. However, as 
Damour and Polyakov speculate [l^, since the dilaton potential is related to the cosmological 
constant, the yet unknown mechanism accounting for the smallness of the cosmological 
constant may also make the dilaton light. Since there is no good theory of the dilaton 
mass, the line labelled "dilaton" in Figure 1 should be terminated a point determined by 
experiment. 
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C. Yukawa moduli 



The Yukawa couplings of the standard model could also depend on moduli which are 
relatively unaffected by Planck scale physics, but obtain a mass due to (low-scale) super- 
symmetry breaking [lj. A Yukawa modulus </> may be coupled as follows to up-type quarks 
and the Higgs boson H: 



£ = \(4>)q L u R H u + h.c. , (7) 

along with analogous terms for down-type quarks and charged leptons. Here flavor indices 
have been suppressed and we for simplicity assume one modulus per coupling. 
Yukawa terms of the form in eq. (J7|) contribute to an effective potential for 0: 

W) = |ll^ At(0)A(0)F2 + W) 

where again F is the fermion-scalar messenger mass-squared splitting, k is a loop-integral 
factor of order 1, N is the number of messenger multiplets, and Vo describes any additional 
unknown contribution to the potential not due to the operator in eq. (|7|). The coupling can 
be expanded around its minimum (0) ~ M, where M is of order string scale, 

and a lower bound on the modulus mass can be obtained from the known term in 



,2 



16^ (A(1)2 + A(O)A(2)) F^ + ^ |0=(0) ^ 



3(16vr 2 ) 3V ' M 2 

For a particular flavor modulus, the coupling = 2m q / \/2vH g where v is the Higgs 
vacuum expectation value of 246 GeV, and H q equals sin/3 for up-type quarks and cos/3 
for down-type, where tan (3 is the ratio of the Higgs' vacuum expectation values. Using this 
result, the Compton wavelength becomes 

F 5 x 10 17 GeV 



= 1050/im x a s 1 



GeV H, 

x ( ^T )( x 



(lOOTeV) 2 M 

,A« 2 A( 2 ) 



A(°) 2 A(°) 



-1/2 



We arrive at an expression larger than reported previously in reference [l( due to a corrected 
numerical factor. 
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The long-range force potential can be determined by relating the scalar 0-quark coupling 
of eq. (JTj) to the scalar 0-nucleon coupling. The fields <p corresponding to up and down 
quarks can generally have different couplings to the proton and the neutron, leading to small 
violations of the equivalence principle. The scalar coupling of the field <fi to the nucleon N 
is again expressed in terms of eq. (j2J) where 

A« M PL (N\m g qq\N) 

The results for the Compton wavelengths and the strengths of the moduli forces relative to 
gravity are plotted in Figure^ The areas are obtained by taking kN = 1, tan (5 = 1, \^ = 0, 
by varying \[F between 30 and 2000 TeV, and by varying A (0) /A {1) x M/(5 x 10 17 GeV) 
between 10~ 2 and 10 2 . 



III. FORCES FROM PARTICLES IN LARGE EXTRA DIMENSIONS 

In theories with large extra spatial dimensions the fundamental scale M* and the observed 
four-dimensional Planck scale M4 = 2.43 x 10 18 GeV are related by 

Ml = m; +2 k 

where n is the number of extra dimensions and V n is their volume. In this framework, the 
standard model gauge and matter content is confined to a 4-dimensional submanifold, and 
the graviton can propagate in all 4 + n dimensions. The scenario provides an alternative 
solution to the gauge hierarchy problem jjj U, 0], as the fundamental scale can be of order 
TeV. Such a paradigm also predicts a modification to Newton's law of gravitation at distances 
nearby and below the length-scale of compactification. 



A. Gravitons in the Bulk 



The modification of Newton's law of gravitation in theories with large extra spatial di- 
mensions has been studied in some detail jiol . l^ . For equal-size extra dimensions and 
toroidal compactification, the volume satisfies V n = (2irR) n = L n , and we have [J| 

R n = 2 x lO^-^mm x (l^X)^ 2 /" 



S 



For the case of two equal extra dimensions (n=2), the radius of compactification is of order 
1-millimeter for M± +n ~ TeV. However for this case, astrophysical bounds require the funda- 
mental scale M 4+ri to be pushed above 1600 TeV [21]. At this scale the two equal-size radii 
are only 2.4 A. Requiring such a high energy scale in turn necessitates more fine tuning for 
the framework to address the hierarchy problem. For the case of 3 extra dimensions M 4+n 
must exceed 60 TeV, and the radius becomes ~ .05 A for equal-size dimensions. These 
limits are derived from Kaluza-Klein gravitons that would be gravitationally trapped and 
remain as a halo surrounding neutron stars j^]. The constraints come from neutron star 
heating via Kaluza-Klein graviton decays. Somewhat weaker limits are also obtained from 
EGRET gamma-ray flux measurements of nearby supernovae and neutron stars. Upcoming 
measurements planned with the GLAST satellite may improve bounds further or lead to a 
new discovery. Both sets of constraints are weakened if the dimensions are of unequal size, if 
there are additional fast decay channels such as other branes for the KK gravitons to decay 
into, or if graviton emission is suppressed as in reference [23j. 

From the four-dimensional point of view, the higher dimensional graviton with momen- 
tum in the extra dimensions appears as a massive particle, leading to a sum of Yukawa 
potentials from the tower of Kaluza-Klein (KK) modes in addition to the massless graviton 
potential. At distances or order R n , only the lowest massive mode contributes significantly 
while the higher modes are exponentially suppressed. For distances r « R n , many modes 
contribute and change the power law dependence of the force from the Newtonian 1/r 2 to 
l/r 2+n . Corrections to the newtonian potential between two masses mi and mi are typically 
parameterized according to the form 

V(r) = -G N ^(l + ae-^) 
r 

where a and A characterize the strength relative to gravity and range of the new force, 
respectively. For distances of order R n or greater, the range A is the inverse of the lightest KK 
mass and the strength a equals its degeneracy [3, Q • At shorter distances, more massive 
KK modes contribute until eventually the power law behavior of the force changes. We 
adopt the convention of references and consider only the leading term corresponding 

to the lightest massive KK modes. For example in the cases of toriodal and spherical 
compactification [3| 

nrU™^^(l + 2n„e^») 
9 



V» n _ sphere = -^1^(1 + {n + l)e i~V-nr, R)) 

for no equal radii of size Rq- 

Due to the stringent astrophysical constraints on two- or three- equal sized large extra 
dimensions, and their even smaller size for n > 3, it is unlikely in this case that the KK 
gravitons can be observed at table-top experiments. However, we stress that these con- 
straints strictly apply to the case of equal extra dimensions. If the extra dimensions are not 
of equal size, it is possible some of the dimensions may be large enough to be detected. To 
illustrate this we plot the cases in (a, A) space for toroidal compactification with one or two 
large radii (amongst possibly many smaller radii). 



B. Gauged Baryon Number in the Bulk 

If in addition to the graviton there are bulk gauge particles, their effective four- 
dimensional gauge coupling g\ can be many orders of magnitude stronger than gravity 
P). As a particular case we consider gauged baryon number B, with the gauge symmetry 
spontaneously broken only on a different submanifold than our own. As discussed in refer- 
ences j^, 24], such a situation can lead to an enormous suppression of the proton decay rate. 
In the following we systematically explore the parameter space for the force strength and 
range. Several details are deferred to the Appendix. 

Since ordinary matter is primarily composed of baryons, the mass of a macroscopic object 
is roughly in proportion to the number of baryons, apart from small effects due to binding 
energy and the electron mass. The expected ratio of the gauge to gravitational forces will 
take the form 

P „2 i 
= r 9auge _ ih_ 1 , x 

9 F grav 4nG N mj 1 ) 

where m p is the mass of the proton. The effective four- dimensional gauge coupling is related 
to the massive 4 + n-dimensional coupling by the volume of the extra dimensions 

1 V d 



dl 9(4+n) 



We can express the (4+n) dimensional coupling in terms of an ultraviolet cut-off scale 
for the gauge theory A which we expect to be of order M*: g\ +n = A~ n 4+n p, where 
0d = ^(d-i)/(27r) d is the 1-loop suppression factor. If the coefficient p is 0(1), this signifies 
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strong coupling in the 4+n dimensional theory, as loop effects become comparable to tree- 
level. Expressing the 4-dimensional coupling in eq. © in terms of A ~ M*, we find the 
baryon-number force can easily reach strengths of 10 6 — 10 8 times gravity, and even higher 
magnitude for strong-coupling with a large number of extra dimensions. 

To avoid conflict with experiment, the baryon number gauge field must acquire a mass. 
If the gauge symmetry is spontaneously broken by a scalar field \ obtaining a vacuum 
expectation value (%), the resulting mass of a gauge particle becomes tra = 9±{x)- Here 
again we assume that \ condenses on a brane other than our own. It was shown in reference 
|2fi| that forces mediated by bulk gauge fields can be exponentially weaker than gravity if 
the bulk gauge symmetry is spontaneously broken on our brane. For (x) = (3 M* with (3 of 
(9(1), the Compton wavelengths are in an interesting range for sub-millimeter experiments. 
A range of predicted parameter space is provided along with a more detailed analysis in the 
Appendix. A portion of the allowed phase space appears in Figures 1 and 2 for comparison 
with other sub-millimeter forces and experimental bounds. 

Astrophysical bounds similar to those discussed for gravitons [i| , [2]]] , can also apply 
for gauge particles in the bulk. In Supernova 1987A, about 10 53 ergs of gravitational binding 
energy was released in a few seconds. One requirement is therefore that the total luminosity 
of Kaluza-Klein particles does not exceed ~ 10 53 erg s _1 . The temperature of the supernova 
is approximately T ~ 30 MeV, and most KK particles are produced with an energy of order 
~ 100 MeV. The constraints on the total luminosity for KK gravitons imply M* > 30 TeV 
for n = 2 extra dimensions [5]. For the case of gauge particles in the bulk, we expect the 
same amount of energy that would have gone into KK gravitons to now produce KK gauge 
bosons, producing roughly the same number of particles. However the rate of production 
for gravitons goes like T n /M 2+n and for gauge bosons like T n ~ 2 /M™, which is more rapid. 
Therefore the constraints on the fundamental scale due to graviton emission for n extra 
dimensions apply for gau ge particle emission with n + 2 dimensions. A similar situation is 
discussed in reference |24| for the case of bulk scalars. This implies the more stringent limit 
of M* >~ 30 TeV for the case of n = 4 extra dimensions. 

Neutron Star Limits. Kaluza-Klein excitations of the gauge bosons can also become 
gravitationally trapped and remain for some time in a cloud surrounding neutron stars. 
Their subsequent decays into photons can in certain cases produce observable gamma-ray 
signals detectable by EGRET. The decay width for such particles can be roughly computed 
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as 

Mi T 

r = iir- < 10 » 

where T is the temperature, typically of order 30 MeV for a supernova. The lifetime becomes 
~ 10 — 10000 years for T from 1 GeV-1 MeV. (The decay width for gravitons goes as 
T ~ T 3 /Mp 1; leading to lifetimes of order 6 x 10 9 years.) However, the situation can be quite 
different depending on the symmetry that is gauged. For B — (X, where nonzero ( denotes 
an admixture of lepton number, there are decay channels into neutrinos. In this case eq. 
(fTU|) is a good approximation and the decays can occur on a time-scale of 10 2 years. Such a 
situation provides little direct observable gamma signal for EGRET. Also KK annihilation 
into positrons is possible, but the resulting gamma rays from positronium annihilation are 
of too low energy to be a useful EGRET source. In the case of pure B, where ( = 0, the 
lifetime given in eq. ()10j) has to be amended since the decay into photons occurs only at 
higher order. The decay width in this case is multiplied by an additional factor of (a/2n) 2 
due to a virtual fermion loop. This increases the lifetime by a factor of ~ 10 6 , making it 
more comparable to the graviton case at ~ 1.3 x 10 s years. The resulting improvement in 
the bound on Jkk as compared with the graviton case is about 50 times. For bulk gauge 
particles, the bound on the compactification scale varies with f KK as M* min a (/kk) 1 ^- 
In this situation, the bounds that applied for n extra dimensions in the graviton case now 
apply to n + 2 extra dimensions. However, even in the case of pure B further limits cannot 
be derived from these neutron star gamma rays and excess heat due to direct re-absorption 
of the KK gauge particles through inverse bremstrahlung, which occurs on a rapid time 
scale and is not loop-suppressed. In this way, there would not be enough remaining KK 
gauge particles to contribute to heating or gamma ray limits. KK re-absorption for the case 
of gravitons was taken into account in reference |22j and in this case the results did not 
appreciably change the graviton limits quoted in reference j^|. As was the case for graviton 
decays, the constraints may be even further weakened if the dimensions are of unequal size, 
or if there are additional fast decay channels such as "photons" on other branes for the KK 
particles to decay into. 
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C. Yukawa Messengers 



An additional possibility is that scalar particles may inhabit the bulk of extra dimensions 



and mediate macroscopic forces, such as the Yukawa messengers considered in reference 24] . 
The messenger fields could be responsible for communicating flavor symmetry breaking from 
other branes to our brane, and for example can attribute the weakness of light generation 
Yukawa couplings to geometrical power-law suppression or exponential suppression due to 
the mass of the messenger fields. In this framework, the vast variation in strengths of the 
Yukawa couplings is recast as a variation in distances to other branes. However, even if the 
messenger fields do not condense on our brane, they can still mediate forces much stronger 
than gravity. If the messenger fields acquire mass due to supersymmetry breaking on our 
wall, their Compton wavelength can be in the sub-millimeter range. As was the case with 
gauged baryon number, the coupling strength for such forces due to the zero mode can be 
large even if the extra dimensions are small enough to make the Kaluza-Klein modes too 
heavy to be detected in sub-millimeter experiments. The coupling strength p ~ v/Mp\ where 
here v is the Higgs vacuum expectation value of 174 GeV and here Mpj is the reduced Planck 
scale of 2.43 x 10 18 GeV. Comparing to gravity, p 2 / (G iv m nucieon) ~ 10 6 . As with the vectors 
described in the previous section, similar astrophysical constraints apply to the scalars in 



the bulk as discussed in reference 



241- 



IV. USING BOSE-EINSTEIN CONDENSED ATOMS TO PROBE (SUB)- 
MICRON DISTANCES 

In recent years the field of atomic interferometry has produced a series of amazing mea- 
surements, including extremely high precision measurements of the acceleration due to the 
Earth's gravity at the level of a part per billion or better. These advances have been made 
possible due to the remarkable techniques developed for trapping and cooling alkali atoms 
(see, for example, the 1998 Nobel lectures of Chu, Cohen- Tannoudji, and Phillips 2jJ). 
Among recent experiments have been a series of atomic-fountain type measurements where 
an atomic beam is launched upwards and allowed to accumulate a phase shift in the Earth's 
gravitational field in a Mach-Zehnder-type interferometer configuration [28]. Similar exper- 
iments have also been carried out to perform sensitive measurements of gravity gradients 
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29[. Such precision techniques can also be used to measure large deviations from Newton's 
constant at short distances. Among the challenges of applying such systems to study grav- 
ity from nearby macroscopic objects is obtaining the optimized beam size and divergence 
necessary to allow a short range interaction to be carried out in a systematic way over a 
long enough time period. Typical spatial extents and velocity spreads are of order mm and 
cm/s respectively, making sub-micron experiments difficult. Experiments involving atoms 
trapped at fixed separation from a source mass surface are in this respect preferable. 

Since its first experimental realization in 1995, Bose-Einstein condensation of alkali gases 
has become a widely growing area of research (see, for example, the 2001 Nobel lectures 
of Cornell, Ketterle, and Wieman 30] or reference |3l| for a recent review). Interference of 
atomic de Broglie waves which develop a relative phase shift due to the Earth's gravity has 
been observed in vertical arrays of trapped Bose-Einstein condensed atoms [32]. The traps 
were located at the antinodes of a laser standing wave, with the trap well depth determined 
by the laser intensity. 

In the following we describe a setup involving arrays of Bose-Einstein condensed atoms 
trapped nearby a surface at the nodes or antinodes of a laser standing wave. The traps are 
thus loaded with atoms in coherent superpositions of states localized at differing distances 
from the surface. For each potential well, the de Broglie phase of the center of mass wave 
function of the atoms evolves according to the interaction potential of the local environment. 
The wave function of atoms localized at different potential wells of the laser can accumulate 
a differential phase shift due to the distance-dependence of the interaction potential with 
the wall. The potential will generally be a superposition of the Casimir-Van der Waals 
potential along with other backgrounds and possibly new short-range interactions. By ad- 
equately subtracting out the Casimir interaction and other background interactions as we 
discuss below, significant improvements can be made over previous searches for new forces 
below 1 micron. In particular the improvements could be 6 to 9 orders of magnitude at 1 
micron, allowing forces of 1000 to 10 6 times gravity to be detected at these distances. Such 
short length scales have been relatively inaccessible to tabletop torsion balance and micro- 
cantilever experiments due to the necessity of having nearby moving macroscopic mechanical 
parts and the unfavorable scaling of the gravitational force with their size. 
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FIG. 3: Proposed experimental arrangement. The standing wave cavity formed by reflection from 
the Au Casimir shield is used to trap BEC atoms in the two potential wells nearest to the periodic 
source mass array. 

A. Experimental Setup and Geometry 

We consider 87 Rubidium atoms prepared in coherent superpositions of states localized 
at the nodes or antinodes of a standing wave of an infrared laser of wavelength Ai as . For 
a laser wavelength that is red de-tuned from the dominant Rubidium D2 line, the atoms 
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are attracted to regions of high laser intensity, corresponding to trapping at the anti-nodes. 
On the other hand, for blue de-tuned light the potential minima occur at the nodes of the 
standing wave. In the first case the potential wells will be centered at distances (Ai as /4, 
3Aias/4,...) from the surface, in the case of normal laser incidence. For blue de-tuning the 
wells are located at distances (Ai as /2,Ai as ,...) from the surface. The effective well separa- 
tion and surface separation can readily be made larger in the case of oblique incidence at 
angle 8, where k\ < k 2 and the effective wavelength determining the trap spacing becomes 
A = Ai as / cos#. Thus a variety of surface and well separations are attainable depending on 
the trap geometry. The atomic well depths can be adjusted by varying the laser intensity to 
overcome the atom-surface interaction potential and the Earth's gravitational field. Trans- 
verse confinement can be achieved for example through the gaussian envelope of the beam 
in the case of red de-tuning, or with additional laser beams. 

For defmiteness we consider a standing wave with an effective trap spacing of A/2 = 
420 nm. We also take the separation of the first trap and the surface to be 420 nm. The 
proposed geometry is illustrated in Figure El The laser is reflected from a 420 nm thick 
shield of gold. The periodic source mass consists of alternating regions of more and less 
dense material, for example gold and silver. Although less dense materials are preferable for 
the Yukawa force contrast, silver is chosen for the similarity of its diamagnetic response to 
that of gold to ameliorate possible problems with magnetic backgrounds. The source masses 
are taken to be 100 microns wide by 100 microns deep by 10 microns tall. The source 
mass can be moved as a whole laterally by a piezo-electric device over several hundreds of 
microns. To reduce temperature changes in the Casimir shield as the source mass moves, 
a very small space is left open directly behind the shield. We consider a population, by 
~ 10 6 atoms, of the first two potential wells closest to the surface, denoted by A and B, 
respectively. The atoms will be spread out in a pancake-like configuration with a transverse 
extent of a few microns. The wells must be loaded with a fixed initial relative phase between 
the two parts of the condensate. The optical wells can be loaded for example by ramping 
up the laser field after evaporative cooling. The population of only the first two wells can 
be achieved perhaps through magnetic or optical techniques to translate the lattice closer 
to the shield wall until the destruction of any extra occupied wells occurs. The de Broglie 
phase of the wave function localized at the first well (A) will evolve more rapidly than that 
of the wave function at the second (B) due to the atom-surface interaction potential. After 
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an interrogation time of l-10s, the laser intensity can be rapidly turned down, allowing the 
wave-packets to escape, spread out, and overlap spatially. The resulting interference pattern 
can then be detected using optical fluorescence. If the accumulated phase difference is due 
to a Yukawa-type interaction with the surface, the phase difference will change depending on 
whether a silver or gold section of the sense mass is positioned behind the Casimir shield. A 
series of such experiments can be performed such that between each experiment the source 
mass pattern is moved laterally behind the screen. Depending on whether a more or less 
dense region of material is behind the shield, the resulting phase shift will display a periodic 
behavior. On the other hand, the Casimir-Polder interaction will be largely the same due 
to the gold Casimir shield. Also, any patch fields that contribute to the interaction can be 
rejected as a common mode. The magnitude of the Casimir-Polder interaction at 420 nm is 
some 8 times the gravitational interaction with the Earth. Due to its gigantic magnitude, 
it is crucial then to consider the finite-thickness corrections to the Casimir potential to 
estimate the potential reach of the experimental setup. 

B. Sensitivity 

Scaling. We compute the difference in the frequency at which the de Broglie phase 
evolves for adjacent populated wells of the laser near the surface due to a Yukawa-type 
potential of strength a and range A. Taking the thickness of the source masses to be r and 2r 
and the atom-wall separation to be r, the potential difference can be roughly approximated 

as 



where we consider a Yukawa potential of range A ~ r and take Sr ~ r, so that we have 



Here m denotes the mass of an individual atom. For example, taking Rubidium-87 and a 
gold wall, the corresponding frequency shift for A = l//m is 1.5a x 1CT 10 Hz. To obtain a 
more precise estimate, we numerically integrate the Newtonian plus Yukawa potential. 

With ~ 10 6 atoms in the condensate, we estimate the minimal detectable phase shift per 
shot as 1CT 3 radians. For a conservative estimate of 1 s interrogation time, the minimal 
resolvable frequency shift is 1.6 x 10 -4 Hz. This corresponds to an acceleration sensitivity of 




(11) 



5V(X) ~ 2aG N mp\ 2 . 



(12) 
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roughly 10~ 7 g, where g is the acceleration of the Earth's gravitational field at the surface, 
where we again take a trap spacing of A/2 = 420 nm. Ultimately, an improvement could be 
obtained by an averaging over 10 4 shots, allowing the phase shift to be detected at the 1CT 5 
level. This along with increasing the interrogation time to 10 s yields a minimal detectable 
frequency shift of 1.6 x 10~ 7 Hz, corresponding to an acceleration sensitivity of order 10 -10 g. 
Obtaining larger interrogation times may be difficult due to loss of coherence from collisions 
or laser instability 3^|. In the following we assume a sensitivity of 1.6 x 10~ 4 Hz. Also, the 
systematic effects we consider in the following sections are not generally problematic at this 
level, however may be more challenging for a measurement with 10~ 7 Hz sensitivity as will 
be discussed. 

Using this minimal detectable frequency shift and the numerical results for the Yukawa 
potential, we generate a plot of the alpha-lambda reach in Figures 1 and 2. Such an experi- 
ment is particularly favorable in the sub-micron length scales, where macroscopic cantilever 
and torsion balance experiments become increasingly more challenging. 

For the geometry described above with a BEC-surface separation of 0.42 microns and 
well separation of 0.42 microns, we evaluate the frequency shift of the potential wells due 
to a Yukawa potential of strength a. For lambda of 1/im, we find 1.2 x 10 _10 a Hz for the 
frequency shift, allowing alpha of ~ 10 6 to be probed. A plot of the projected sensitivity on 
alpha-lambda space appears in Figures 1 and 2 as a solid green line. We also indicate the 
projected improvement possible with the sensitivity taken to the 10~ 7 Hz level as a solid 
red line. In Figures 1 and 2, the green curve is calculated for a surface separation, trap 
spacing, and shield thickness of 420 nm with 1 second interrogation time and a single shot. 
In Figure 2 the red curve is also computed at a well separation, surface separation, and shield 
thickness of 420 nm. The red line in Figure 1 has been computed for a shield thickness, trap 
spacing, and surface separation scaled from 420 nm to 600 nm. Here we assume 10 seconds 
of interrogation time and average over 10 4 shots, corresponding to 10~ 7 Hz sensitivity. The 
larger scale of 600 nm becomes advantageous for measurements near A = 1/im as it can 
better suppress the Casimir background. We note that for these geometries the sensitivity 
levels off above 1 times gravity at large lambda. By scaling the geometry and trap spacing 
together, the estimated sensitivity roughly follows Eq. (fT2"j) so that lx gravity is achievable 
at ~ 30 microns. These larger length scales are also accessible by upcoming micro- cantilever 
or torsion oscillator experiments. Therefore in figures Q and E] we have emphasized the reach 
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at the micron scale and below. 



C. Systematics 



1. Atomic interactions 



Atomic interactions in general will also produce differences in the chemical potential 
of the two parts of the condensate, leading to relative phase differences that would occur 
in addition to those caused by the differential gravitational potential. As one possible 
solution, care could be taken to have an approximately equal number of atoms in the two 
clouds, perhaps by using 2-d lattice co nfigu rations, where arrays of single atoms are confined 



transversely as well as longitudinally 33]. Alternatively the Feshbach resonance could in 
principle be used to highly suppress atomic interations j34J] . For example, it has been shown 
that the s-wave scattering length can in this way be tuned over several orders of magnitude 
and set effectively to zero, thus turning off the chemical potential due to atomic interactions. 



2. Thermal fluctuations 

Finite temperature fluctuations of the BEC phase may present an additional experimental 
challenge, common to interferometry using BEC. For example, such effects have been studied 
theoretically and experimentally 35j] for highly-elongated condensates where they are shown 
to become problematic, as the system becomes quasi-one-dimensional. Care should be taken 
to keep a three dimensional nature to the condensate to minimize this effect. Also a finite 
temperature in the surrounding materials may pose an experimental challenge. Thermal 
currents in metals in the local environment can produce magnetic field fluctuations which 

n 

can limit the condensate lifetime [36J. Such fluctuations can be minimized by lowering the 
system temperature, using lower-conductivity materials, and minimizing the thickness and 
transverse area of the reflecting shield. 

3. Casimir Background 

Finite Thickness and Conductivity. The frequency shift due to the bare Casimir 
force is quite large (on the order of 8 kHz). However, in the case of infinite conductivity, the 
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differential frequency shift between regions of metal of varying thickness is zero. In practice, 
the finite conductivity of the metal as well as its finite thickness has to be taken into account. 
It was shown in reference j3] for the case of two metal walls that the differential Casimir 
background due to differences in thickness rapidly dominates the Yukawa force as thickness 
and plate separation decreases, thus making measurements of gravitational strength Yukawa 
forces difficult below a few microns. The geometry studied consisted of a semi-infinite probe 
mass of finite conductivity separated by a distance D from a finite conductivity source mass 
of thickness D. The result for the Casimir force and a gravitational strength Yukawa force 
of range D was then compared to the case source mass thickness 2D. It was shown that the 
differential Casimir background due to differences in thickness becomes comparable to the 
differential Yukawa force at distances of about 3 microns and rapidly dominates the Yukawa 
force below this length scale. In this work we show that by replacing the metal probe mass 
with a dielectric, or with an atom in particular, the differential Casimir force is considerably 
smaller, making the domination over the Yukawa potential less severe at sub-micron lengths. 
Following references 3| , [3 we employ a reflection-based model for computing the Casimir 
force between two walls. We obtain similar results to reference B[ for metallic walls. Also, 
we obtain a reduction factor rjp, (defined precisely in reference 2zl) which describes the 
reduction of the Casimir force at small separations and for finite conductivity, that agrees 
well with the results of that reference in the case of semi- infinite walls. The parameter rjp 
is defined as 

F c = VfF p (13) 
where Fp is the perfect conductor result. The integral expression 

120 r° ,„ - r K r 2 n 



dKK2 L dn V^i (14) 

gives rjp in terms of the reflection amplitudes r p which can depend on the polarization, con- 



ductivity, frequency, and wall thickness |37[. Here K = kL and f2 = are the wavenumber 
and frequency measured with respect to the cavity length L. 

The Casimir force between the mirrors can be expressed in terms of the imaginary part 
of the dielectric function of the walls. For metallic walls, the dielectric function is 

e(iu) = 1 + 
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where we assume a Drude model for the metals, and ui p and 7 are the plasma frequency 
and relaxation frequency, respectively. The magnitude of the Casimir force changes rapidly 
as the length scales of thickness and separation approach the plasma wavelength X p . The 
expression for the reduction factor rjp which describes the fraction of the perfect metal 
Casimir result can be written as an integral over all frequencies and wavenumbers. In 
particular, the low-frequency response of the dielectric function of metals diverges as uo — * 0. 
By replacing one of walls with a dielectric e, the low-frequency response becomes weaker 
and contributes less to the Casimir force. As a consistency check, we evaluate the expression 
for a dielectric wall and metal wall using both the reflection model of reference j^J an d by 



381 ] . The two numerical 



numerically integrating the zero-temperature Lifshitz result directly 
calculations agree to a part in 10 5 . After demonstrating the equivalence of the two models 
for the semi-infinite case of dielectric and metal walls, we proceed with the reflection model 
to study the finite thickness dependence. We find the result is less sensitive to the thickness 
of the source metal, which improves the situation considerably for new force detection below 
3 /mi. The reduction factor for the Casimir Force rjp for the two walls is computed for the 
case of wall thickness D and 2D, for wall separation L. We define the quantity 

(2D,L=D) _ (D,L=D) 

= * un^ (15) 

Vf 

which expresses the fractional differential Casimir force for the source walls of different 
thickness. We list a table of values of Arj for probe walls of metal and dielectric materials 
below. We also note that the interaction between two dielectric walls is much more sensitive 
to their thickness, due to the lack of screening present in metals. For example, we find Arj 
can be as large as 1 percent for D = 1 /zm and e = 10, which compares rather poorly with 
the metal-dielectric case. 

To achieve the limit of the atom-wall interaction, we consider rarifying the dielectric 
medium. This technique was used by Lifshitz to derive the individual atomic Van der 
Waals potential. We also add the dominant resonance for Rubidium D2 line. The dielectric 
function satisfies 

e(iu) = 1 + 4:irna(iuj) 

where a is the dynamical polarizability and n is the number-density of atoms. Considering 
the dominant D2 line at 780 nm with an oscillator strength of nearly 1, we approximate the 
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TABLE I: Estimates for the differential Casimir force between a gold source mass of thickness 
D versus 2D, and a semi-infinite probe mass of varying materials. Smaller values indicate less 
sensitivity to the differential thickness. 



D 

(/im) 


^??metal— metal 


^^/metal— dielectric 

e = 100 


^^/metal— dielectric 

e = 1.001 


.1 


9 x 1(T 5 


3 x 10~ 5 


2 x 10~ 5 


.3 


5 x 1(T 6 


1 x 10~ 9 


9 x 10~ n 


.6 


2 x 1(T 6 


1 x 10~ 10 


4 x 10" 12 


1 


8 x KT 7 


1 x 10~ 10 


4 x 10~ 13 



polarizability as 

OiitUJ) - 



The role of the parameter e(0) — 1 is now played by the quantity Annao- For 87 Rb, we have 
ao = 2.7 x 10~ 23 cm 3 , and even for a high number density of order n ~ 4 x 10 16 cm -3 we 
obtain e(0) — 1 of ~ 10 -5 , indicating even more favorable scaling of A?7 than shown in the 
Table for 1(T 3 . 

It remains to estimate the minimal detectable alpha due to a Yukawa potential limited 
by the differential Casimir force. The Casimir potential due to the atom at distance r from 
an infinitely conducting surface can be written 

Uc = -— ^- (16) 
° 8tt r 4 V ; 

In practice we expect corrections due to the finite conductivity and the dynamical polariz- 
ability of the atom. The finite conductivity correction is less than a factor of 2 reduction for 
the length scales of interest. The equation is strictly valid in the limit of large separa- 
tion, at length scales greater than X/2n where now A is the wavelength which contributes to 
the atom's polarizability. For lengths below this scale, the Casimir screening due to retarda- 
tion becomes less effective and the power law changes to 1/r 3 corresponding to the van der 
Waals interaction. For Rb, the dominant wavelength is 780 nm and since we are interested 
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in length scales above 100 nm, the form in Eqn. (JTfij) is a reasonable approximation. We 
now multiply the parameter At] by the difference in Casimir potentials of the wells obtained 
through Eq. (|16|) to obtain an estimate of the difference in potential shift due to the metal 
walls of the two thicknesses. The resulting A5Uc we compare with A^f/yukawa- (Here A sig- 
nifies the change from a region of thickness D to 2D, and 5 signifies the different locations of 
the condensates A and B. In Figure 0] we abbreviate both by d = A5.) In order to illustrate 
the scaling, we consider a set-up with equal length D for A, the well spacing, and the atom- 
wall separation, and compare regions of metal thickness D and 2D. We display the results 
in Figure |U and take a Yukawa force of gravitational stength, a = 1. For comparison we 
include as a dotted line the estimate adapted from reference |7| which illustrates the scaling 
of the differential a = 1 Yukawa and Casimir forces for a metallic probe wall in place of the 
atoms. The situation for the atom-wall setup is more optimistic than that of the metallic 
wall- wall setup (shown in Fig [5] of reference 0]) by several orders of magnitude between 
200 nm and 3 microns. In Figure 0] there are two separately labelled vertical axes shown 
since the atom-wall Yukawa detection limit due to the differential Casimir background is de- 
termined by the ratio of the potential differences, whereas for the metallic wall-wall system 
the relevant quantity is the ratio of the forces. This does not prevent a direct comparison 
of the a- reach of the two systems however. For example, we see from the figure that a of 
10 3 can be reached at .6 microns in the metal-atom case and at about 1.5 microns in the 
metal-metal case. 

Temperature Effect. The temperature dependence of the Casimir force becomes quite 
weak for low temperatures and small separations. However, it is important to estimate the 
temperature-dependent contribution since the full magnitude of the Casimir force is so large. 
A perturbation approach has been developed in reference j39j. We employ the parallel metal 
plate expression valid for separations below 2 microns as a function of temperature T: 

A /m . , 4 2515 3C(3) X . , x 

where t = T e g/T, T e g = hc/2akB, 5 = \ p /2tt, and a is the plate separation. Taking a 
separation of a = 840 nm, a gold surface of X p = 136 nm, we find for a room-temperature 
of 295 K, the parameter t = 4.6. This amounts to a fractional change in the Casimir force 
Ar/p = .0012At and writing At as t^r, we have 

ST 

At? f (295K) = 0.0055 x — ... (18) 
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D (microns) 



FIG. 4: Estimate for the Yukawa detection limit due to Casimir-Van der Waals background for 
Yukawa range A = D and strength relative to gravity a = 1 is shown as the black solid line for 
the atom-wall system. The well spacing and atom-wall separation are also D, and we compare 
regions of metal thickness D and 2D. The metal- metal case studied in reference [jj is shown as a 
red dotted line for comparison. 



which can be a significant effect. If the entire reflecting surface changes temperature by 
10~ 6 K during the measurements, the frequency shift is at the 10~ 7 level. For the proposed 
initial parameters the sensitivity to wall temperature is at the ~ mK level. However, most 
of the Casimir force results from the wall material closest to the atom, which we suggest to 
take as a shield of uniform material, and have the pattern of varying density source masses 
separated from this material by a small distance to prevent heating upon motion. Even 
though a small temperature gradient could be supported across the shield wall, it is unlikely 
that the shield surface temperature will vary periodically with the source mass density, as 
the two are not in direct thermal contact. A periodic temperature gradient in the source 
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mass distribution itself changes the expected Casimir force very little, as its effect comes in 
only at the level of A77 due to the thickness of material calculated in the previous section. 
As an additional handle, the temperature of the surfaces can be controlled externally to 
quantify the effect experimentally. Also, decreasing the temperature improves the situation 
considerably, yielding a sensitivity to a full wall temperature change of ~ 1CT 3 K at 77 K 
and only to ~ 10 K at liquid helium temperature for frequency shifts of 1CT 7 Hz. 

Isotope Effect. A unique feature of atomic systems in contrast to macroscopic objects is 
that the electrical properties and mass of the atom can be toggled in a precise way by taking 
advantage of other stable atomic isotopes. For example, by using 85 Rb, and comparing to 
an experiment done with 87 Rb, one expects the force to change at the 10~ 2 level, while 
the Casimir force changes only at the 10~ 4 — 10~ 5 level. To verify this claim, we compute 
the Casimir force for each isotope according to the Lifshitz model, and simply change the 
wavelength of the dominant transition from 780 nm to 780.1 nm. Although this technique 
decreases the sensitivity to alpha by 2 orders of magnitude, it may be useful for doing 
measurements at around 100 nm separation from the surface, where the differential Casimir 
force rapidly dominates over the Yukawa force due to the finite plasma wavelength of the 
metal. 



4- Magnetic and Other Backgrounds 

Magnetic Susceptibility. Local magnetic field gradients can be caused by the variation 
in magnetic susceptibility of the two source mass materials. For a background magnetic field, 
e.g. from the Earth, the induced magnetic dipole moment in the materials produces a field 
which varies in proximity to the two materials. The induced magnetization in a paramagnetic 
or diamagnetic material satisfies 

M — - Xn \ B 

(1 + Xm)^0 

Now the induced field due to the magnetized materials we roughly approximate using the 
expression for a magnetized sphere 

2fi Ma 3 2 X m D a 3 

B{z) ~ ~ JTTxm) 

— * 

where a is the sphere radius and B is the background magnetic field responsible for the 
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induced magnetization. We find the differential field to be 

A „ . , 2Ay „ a 3 ,Az, 
(1 + x) z 3 z 
2 

~ AyB 

100 A 

where in the last line we assume Az = 420 nm and we take the radius to be a = 50/im. 
Silver and Gold are both diamagnetic, with a differential susceptibility of .6 x 10~ 5 . We 
estimate AB ~ 10~ 7 .B. The frequency shift due to magnetic fields for Rubidium is 1.4 
MHz/Gauss, so to obtain 10~ 4 Hz resolution only requires a magnetic field shielding of 
B < 1 mG. One can further alleviate the constraints on background magnetic fields by 
choosing materials that have more similar magnetic susceptibility either in pure form or 
through selective doping. Although not problematic for the initial proposed parameters, 
extending the Yukawa sensitivity down to the ~ 10~ 7 Hz will require more extensive magnetic 
shielding or precisely tailored alloy or doped materials. 

Gravity as a background. Although heavy nearby objects can be easily detected, this 
is not expected to be problematic since they in general cannot exhibit the periodicity of 
the source mass pattern. We note that in order to avoid acquiring a differential background 
signal at the 10~ 7 Hz level due to the gravitational attraction of the proof masses themselves 
(which is only power-law suppressed and so remains significant at distances much greater 
than the A of interest below 1 /im), it is necessary to limit the vertical extent of the proof 
masses to be less than approximately 100/xm. 

Finally, we list a number of other systematic backgrounds which are not expected to be 
problematic due to common-mode rejection. They include patch field effects on the surface 
of the reflecting metal shield, the roughness of the surface of the shield, the background 
Earth's gravitational field. 

We conclude this section with a table summarizing the expected frequency shifts of se- 
lected systematics as they compare to the Yukawa signal. Improvements beyond the level we 
discuss may be attainable by tailoring materials to have more similar conductivity, magnetic 
susceptibility, and by going to low temperatures. 
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TABLE II: Summary of selected systematic effects. The three backgrounds: Newtonian, A Casimir, 
and magnetic are evaluated for a surface separation, trap spacing, and shield thickness of 420 nm. 
Scaling these distances from 420 nm to 600 nm causes the differential Casimir signal to drop below 
the 10 -7 Hz level. 



Yukawa signal 




A 


A/(Hz) 


1 fim 


1.2 x 10~ 10 a 


.6 ^m 


3.6 x lO -11 a 


.3 ^m 


3.4 x 10~ 12 a 


.1 fim 


1.8 x 10~ 15 a 


Background 




Newton 


6.5 x 10" 9 Hz 


Casimir 


8283 Hz 


Ar) 


7 x 10 -11 


ACasimir 


5.7 x 10~ 7 Hz 


Magnetic 


3 x 10~ 5 Hz x£\B (mG) 



V. DISCUSSION 

The search for short-distance modifications to Newtonian gravity has rapidly expanded 
over the past seven years. We have illustrated several possible examples of sub-millimeter 
physics that may occupy a vast amount of still unexplored parameter space. We have also 
described an experimental technique involving interference from arrays of Bose-Einstein 
condensed atoms that could extend the search by several orders of magnitude below a micron. 
Such techniques, if successful, allow access to an area of phase space that is complementary 
to the super-micron reach of upcoming torsion oscillator and micro-cantilever experiments 
and could lead to exciting discoveries. 
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FIG. 5: Interaction strength and range for gauged Baryon number in the bulk, a = 1 corresponds 
to a force of Newtonian gravitational strength. For illustration, the case of (3 = 1 is shown as the 
green right-most parallelogram and (3 = 2 is shown in yellow. The ranges for strong and weak 
coupling in 3 and 6 extra dimensions are shown. The upper limit shown for the strong-coupling 
region in the case of 3 extra dimensions terminates at smaller alpha as shown. The weak-coupling 
lower boundary is identical for both cases. 

VI. APPENDIX 



In this appendix we estimate the range and strength of forces due to gauged baryon 
number in the bulk. The mass of a macroscopic object is roughly in proportion to the 
number of baryons, and the expected ratio of the gauge to gravitational forces takes the 
form 



a: 



gauge 



grav 



i 



47T G N m 2 v 



(19) 
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where m p is the mass of the proton. As discussed in the text, we take the 4+n dimensional 
gauge coupling as gf +n = A~ n 4+ „p, where A is the ultraviolet cutoff scale and @d = 
fi(d-i)/(27r) ci is the 1-loop suppression factor. If the coefficient p is 0(1), this signifies strong 
coupling at the scale A in the 4+n dimensional theory, as loop effects become comparable 
to tree-level. Using ild-i — 27r ci//2 /r(|), we find the familiar 4-dimensional loop factor 
04 = 167T 2 . For three extra dimensions the factor becomes 07 = 15 • 87r 4 and for six extra 
dimensions W = 12 ■ 10247T 5 . The effective 4-d gauge coupling satisfies 

1 = Vd 

,M 4 2 1 1 ,A d , 



y Mj V n . d Mr d d+ AP K Mf 

where in the second line we have assumed that the gauge bosons propagate in d < n of the 
extra dimensions. For definiteness, we consider the case of gauge bosons living in all extra 
dimensions (n=d), and the expression simplifies to 

M 2 M n 

s 4 2 = (^)0, + 4-p(fr). (20) 

so then the ratio of forces becomes 

<* = T-7±-> < 21 > 

= (^ )25 0^(^). (22) 
v m p 2 ' 4tt v A n ' v ' 

The factor of 25 appears as the square of the approximate ratio between 1 / y/G^ and the 
reduced Planck mass. For M* of order TeV, already the force becomes ~ 1 million times 
gravity. If the gauge symmetry is spontaneously broken by a scalar field x obtaining a 
vacuum expectation value (x), the resulting mass of a gauge particle becomes 

m A = 5-4 (x)- (23) 

Here again we assume that x condenses on a brane other than our own. It was shown in 
reference that forces mediated by bulk gauge fields can be exponentially weaker than 
gravity if the bulk gauge symmetry is spontaneously broken on our brane. 
We parameterize the vacuum expectation value of x as 

< X >=P-A (24) 
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where beta is a numerical coefficient of order 1, so that the Compton wavelength is written 

1 M 4 1 A"/ 2 

If the scale A is somewhat smaller than M*, the Compton wavelength becomes shorter and 
the strength a g increases. For example, if one imagines new physics arising from string 
theory at scale M s , we have the relation 

M 2+n 

so that 

M s = M*(^) 1/(2+n) < M*. (26) 
where we have identified the Lagrangian 

/ rf 4 xv^^M 2+n K = / d^xyf^gllMl (27) 

with a Type I string Lagrangian j^J 



/ 



<? 2 



If A is associated with physics at the string scale, it is possible the coefficient (3 is greater 
than 1. 

In order to illustrate the phase space encompassed by these forces, in Fig. 1 we plot a 
versus A^, where for simplicity we set A = M* and we vary the parameter p from 1/(10 x 
0^+4) for weak coupling, to 1 for strong coupling. We vary M* from 1 TeV to 100 TeV, and 
show f3 between 1 and 2. For smaller (3, the predicted region can be extended horizontally 
to the right until conflicting with experimental observation. 

Finally we note the gauge force we discuss is due strictly to the zero-mode, and therefore 
its range is not strongly limited by the size of the extra dimensions. In the case that any of 
the extra dimensions has a large enough compactification radius, the lightest KK modes of 
the gauge particles may also make a contribution, though we do not include this explicitly 
in Figures 1 and 2. 
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